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Abstract 

Using the ATLAS 4.7 fb" 1 data on new physics search in the jets + tf)r channel, 
we obtain new limits on the lighter top squark (t\) mass assuming that it is the next 
to lightest supersymmetric particle (NLSP). If the decay t\ —> cxi dominates and the 
production of dark matter relic density is due to NLSP - LSP co-annihilation then the 
lower limit on m~ is 240 GeV. The limit changes to 200 GeV if the decay t\ — > bWxi 

t> 1 

dominates. Combining these results it follows that t± NLSP induced baryogenesis is 
constrained more tightly. 
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Models with Supersymmetry (SUSY) [1] are very well motivated extensions of the Stan- 
dard model (SM) of particle physics. The LHC experiments [2, 3] with 4.7 fb -1 data have 
severely constrained the masses of the strongly interacting supersymmetric particles (spar- 
ticles) - the squarks and the gluinos [1]. However, a subset of these sparticles may still be 
rather light. An especially interesting possibility is the light stop scenario. In this case only 
one of the squarks - the lighter superpartner of the top quark or stop (ii) - is light and 
within the reach of the ongoing LHC experiments along with the electroweak sparticles (the 
sleptons and the electroweak gauginos) 3 . 

The light stop scenario is theoretically well-motivated. Due to a large mixing term - 
driven by the top quark mass - in the stop mass matrix, one of the mass eigenvalues (i.e 
m^) may turn out to be much smaller than the masses all other squarks. This result is 
model independent. Additional suppression of may arise due to the underlying SUSY 
breaking mechanism. In the minimal supergravity (mSUGRA) model [8], e.g., all squarks 
have a common mass at the Grand Unification Theory (GUT) scale (Ma)- Yet at the 
weak scale may be suppressed relative to the masses of the other squarks due to renormal- 
ization group running between Mq and the weak scale driven by the large top quark Yukawa 
coupling. In this paper,however, we consider the light stop scenario at the weak scale as 
a phenomenological model without invoking any specific mechanism for the suppression of 
mj . Light stop squarks are cosmologically interesting since they can provide an explanation 
of baryogenesis [9]. This happens if 120 ~ < m t . 

If t\ is the next to lightest super particle (NLSP) then another cosmologically interesting 
possibility opens up. The observed dark matter (DM) [10, 11] relic density of the universe 
[12] may be produced by the co-annihilation of t\ with the lightest supersymmetric particle 
(LSP), the DM candidate [13]. We shall consider models where the lightest neutralino (xi) 
is the LSP. For this co-annihilation, however, Am = mr — m~o has to be rather small. The 
novel idea that the light stop scenario plays the major role in both baryogenesis and DM 
relic density production has also been emphasized in the literature [14]. In this work we 
shall focus on the possible constraints on the stop NLSP scenario for both large and small 
Am from the published LHC 7-TeV data [2] corresponding to integrated luminosity C = 4.7 

3 It has recently been emphasized [4, 5, 6] that if the squarks and the gluinos are compatible with the 
LHC mass bounds, the data is consistent with a wide variety of relatively light electroweak sparticles and 
they may have masses just above the corresponding lower limits from LEP [7]. 
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fb _1 . Similar analyses were done for £ = 1 fb _1 data in [15] 

Before starting the main discussion it is worthwhile to review the existing bounds on the 
stop NLSP from the LEP and the Tevatron. The Stop NLSP has the following decay modes: 
(i) ti -> cxi (the loop induced mode), (ii) ti -)• bWxl, (hi) h -)• tx\ and (iv) t x bff'Xi 

The Branching ratios of the modes (i) and (iv) are negligible if either (ii) or (iii) is 
kinematically accessible. The modes (i) and (iv) can in principle compete with one another 
[16], but the latter is negligible over most of the parameter space if tan (3 is moderately large 
(> 7) [16, 17] and will be ignored in this paper. If the t\ — Xi co-annihilation occurs then 
only decay mode (i) is allowed. 

Most of the existing bounds on stop NLSP are derived by assuming that the loop induced 
decay (i) occur with 100%. The expected signal is (at least) 2 jets + fir. For m^r within the 
kinematical reach of LEP, the bound is sensitive to small values of Am as well and is given 
by mq > 100 (95.2) GeV for Am > 20 (5) GeV for no mixing (9f — 0) [7]. The numerical 
change in the limit is not drastic for other values of % On the other hand the Tevatron 
experiments are sensitive to larger m-jr but are valid for rather larger Am. The best bound 
is m^ > 180 GeV for m~o = 90 GeV. Thus stop induced baryogenesis is disfavoured for Am 
~ 90 GeV. For smaller Am the efficiency of the kinematical selection requiring high jets 
becomes rather poor due to the modest energy release in the decay process and the limit 
becomes weaker. The possibility of stop NLSP induced baryogenesis, therefore, remains open 
for small Am. It is, also, fair to conclude that the possibility of t\ - LSP co-annihilation is 
only excluded for m^- ~ 100 from the LEP bound. 

The ATLAS or the CMS collaboration has not yet published any dedicated search for 
the pair production of the stop NLSP. However, this pair production followed by the decay of 
each stop in channels i) - iv) would certainly contribute to the jets + missing signature 
carefully looked into by the LHC collaborations. The non-observation of this signal can 
potentially constrain the light stop scenario. It has, however, already been noted that the 
£ « lfb -1 data was sensitive to the stop sector [15]. However, we shall show below that the 
limits are much stronger for the £ = 4.7 fb _1 data. This large £ accumulated at the LHC 
enables one to constrain the stop NLSP scenario meaningfully for smaller Am and directly 
address the issue of Stop NLSP - LSP co-annihilation as a DM producing mechanism at a 
hadron collider. 

We have computed the sparticle spectra and decay branching ratios (BRs) in the minimal 
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supersymmetric extension of the standard model (MSSM) without any assumption regard- 
ing the soft breaking parameters using SUSPECT (v2.41) [18] and SDECAY [19]. Using 
PYTHIA (v6.424)[20] Monte Carlo (MC) event generator we have generated the jets + fa 
signal stemming from t\t\ events only. The next to leading order (NLO) cross-section for the 
htl pair production have been computed by PROSPINO 2.1 [21] with CTEQ6.6M PDF [22]. 
The normalization and factorization scale is chosen to be vnq . We have taken the pole mass 
of the top quark (running bottom quark mass evaluated in the MS scheme) as m t (nib) = 
173.2 (4.25). We have computed the DM relic density using micrOMEGAs (v.2.4.1) [23]. 
The observed DM relic density (Qh 2 ) in the universe measured by the WMAP collaboration 
[12] is given by Vth 2 = 0.1126 ± 0.0036. If 10% theoretical uncertainty is added [24] then 
the DM relic density is bounded by 0.09 < Vth 2 < 0.13 at 2a level. 




Figure 1: The distribution of Nj e t (in the left panel) and pr (in the right panel) of the hardest jet for the 
Benchmark points BP1 and BP2 (see text for details). 

ATLAS Collaboration has updated their results for SUSY searches in jets + fa + 0/ final 
state for C = 4.7 fb _1 [2]. They have defined five inclusive analysis channels labelled as A to 
E according to jet multiplicity (Nj = 2 to 6). Further, the Nj = 2 sample is divided into two 
channels - A and A'. The channels and the details of kinemetical cuts for each of them are 
given in Table 1 of [2]. Depending upon the final cuts on the observable m e //(incl.) (defined 
as the scalar sum of all jets with > 40 GeV and fa) each channel is further classified as 
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'Tight', 'Medium' and 'Loose'. In this way they have finally presented the results for 11 signal 
regions and have constrained any new physics model in terms of upper limit on the number 
of events Nbsm or effective cross-section (Tbsm/R)- The observed upper limits on Nbsm at 
95 % Confidence Level (CL) for signal regions SRA-Tight, SRA-Medium, SRA'-Medium, 
SRB-Tight, SRC-Tight, SRC-Medium, SRC-Loose, SRD-Tight, SRE-Tight, SRE-Medium, 
SRE-Loose are 2.9, 25, 29, 3.1, 16, 18, 58, 10, 12, 12, 84 respectively [2]. 

We have adopted the selection criteria for different signal regions used by the ATLAS 
Collaboration [2]. Then we have checked whether the number of events from the tit\ pair 
production followed by the decays (i) - (iii) of both the squarks exceed the corresponding 
upper bound for at least one of the signal regions to obtain the exclusion contour in the vnq 
- m~o plane. 

If the BR of the decay mode ti — > cXi is 100 % then the signal will be 2j + fsp. In this 
case SRA-Medium or SRA'-Medium, specially designed for compressed spectra [2] is very 
effective in excluding the regions of the — m~o space. For small Am, SRA'-Medium is 
the most sensitive one. On the otherhand, high jet multiplicity regions chosen by ATLAS 
are optimised for long decay chains [2]. For t\ — > bW\\ or tXi decay, events will contain 
at least 2, 4 or 6 jets. Due to the lepton veto 2j signal will not be sensitive at all. In our 
analysis we have also checked that when BR of t\ — > tx® is 100 % then the exclusion mainly 
comes from SRE-Medium (Nj = 6) channel. Again for t\ — > bWxi effective signal regions 
are either SRC- Medium or SRE-Medium. The distributions of the number of jets (Nj) and 
the pt of the hardest jet in the signal are shown in Fig. 1 for the two representative point 
BP1 and BP2. For both of them m 7 is 200 GeV. But m~ = 180 and 100 for BP1 and BP2 
respectively. Both the distributions are drawn with nominal cuts and normalized to unity. 

Fig. 2 depicts the regions in the - m~o plane excluded by the ATLAS jets + Ifa data. 
In the region between the mr = m c + m~o and the mr = m^ + my/ + m~o lines (Region-I), 
the loop decay (i) dominates. We have restricted ourselves to (mj^min > m c + m~o + 2 
GeV. The same increment of 2 GeV has been introduced at each point representing a cross 
over from one region to another. In Region-II bounded by the lines = m;, + my/ + m~o 
and = m t + m~o lines, the three body decay channel (ii) dominates. In regions I and 
II the parameter space to the left of the curved line is excluded. In Region-Ill below the 
line = m t + m~o , the two body decay mode (iii) is the only significant decay mode. 
The parameter space bounded by the dashed line and the x-axis is excluded. The following 
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Figure 2: Exclusion limits in the m~ - m~ plane for different decay modes of the t\ - NLSP. The x marked 
points represent the parameter space where t\ — xi co-annihilation is viable. 

points may be noted : 

a) For very small Am the jets in the signal essentially come from initial state and final 
state radiation. This result is, therefore , crucially dependent on the parton showering model 
in Pythia. In contrast for Am = 25 - 30 GeV, which is the case in the co-annihilation region, 
even the parton level pt distribution develops a long tail. Our results for this region is seem 
to be more reliable. In the co-annihilation region vrtqr ~ 240 GeV is disfavoured. In Region-II 
the best limit is mr ~ 230 GeV for m~o ~ 150 GeV. 

b) The t\ NLSP induced baryogenesis requiring 120 ~ mq < m t seems to be disfavoured 
in the co-annihilation region. In Region-II vnq ~ 200 GeV and hence, t\ induced baryogenesis 
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is also disfavoured. In Region-Ill it is forbidden by definition. 

c) Significant constraints are obtained in decay channels (ii) and (iii). The strong con- 
straints obtained even for small Am are due to the large £ accumulated at the current LHC 
7 TeV experiments. Most of the above bounds come from ATLAS selection criteria requiring 
two high px jets in the final state . That is why the constraints in Region-Ill are rather 
modest. Here the final state contains either large number of jets and/or at least one high pt 
lepton which is vetoed. 
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